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Abstract. Following the derivation of the trace formulae in the first paper in this 
series, we establish here a connection between the spectral statistics of random regular 
graphs and the predictions of Random Matrix Theory (RMT). This follows from the 
known Poisson distribution of cycle counts in regular graphs, in the limit that the cycle 
periods are kept constant and the number of vertices increases indefinitely. The result 
is analogous to the so called "diagonal approximation" in Quantum Chaos. We also 
show that by assuming that the spectral correlations are given by RMT to all orders, 
we can compute the leading deviations from the Poisson distribution for cycle counts. 
We provide numerical evidence which supports this conjecture. 



1. Introduction 

The present paper is the second in this series, where we aim to estabhsh a rigorous 

connection between the spectral fluctuations in the spectra of random regular graphs 
and the predictions of Random Matrix Theory (RMT). In the first paper [1] (to be 
referred to as /) we provided the necessary definitions and facts about graphs and we 
shall use the same notations here. Suffice it to say that wc deal with the ensemble Qy^ 
of d-regular graphs on V vertices, and we study the spectrum {fxjjj^-^ of the adjacency 
matrix A, from which we excluded the trivial eigenvalue fxy — d. The spectral density 
is defined as 



In what follows we shall be interested in the large V limit, and in most cases the 
replacement of — 1 by will be justified. We shall do this consistently to simplify 
the notation. 




v-i 



(1) 
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In / we also defined the ensemble of "magnetic" graphs, Qy^^ where the adjacency 
matrix of each member of the ensemble Qv,d is decorated by phases 

^ Aje'^'-^ ; Xj,i^-Xi,j (2) 

and the phases are independent random variables distributed uniformly on the unit 
circle. Here the entire spectrum of M is considered and 

P^''\n)^l^j2^{n-Hj). (3) 

i=i 

In / we prepared the tools needed for our purpose, namely trace formulae. In the 
sequel, we shall summarize the absolutely necessary information about trace formulae 
required to make the present paper self contained. 

1.1. The trace formula - a short reminder 

Trace formulae express the spectral density of the adjacency matrix as a sum of two 
contributions, = ^ ^(/^ — l^j) = {p{f^)) + p(/x) where ( ■ ) stands for the 

ensemble average. Here, For both ensembles, {pip)) is the well known Kesten- McKay 
expression for the mean spectral density [4, 5] 

Pkm(p) = ,hni (p(/x)) = < . (4) 

t for \^\ > 2^fd^ 

Note that the Kcstcn-McKay density depends explicitly on the degree d. p{fi) is the 
fluctuating part of p{p), with {p{p)) = 0. To simplify the notation we omit reference to 
the d dependence. In the limit d^ 1, pkm approaches Wigner's semi-circle distribution 
which characterizes the canonical Gaussian random matrix ensembles. The fluctuating 
parts are expressed as infinite sums over Chebyshev polynomials (of the first kind) Tn{x) 
with coefficients which depend on cycles on the graph. The trace formulae take similar 
forms for the two ensembles, but with different coefficients. 

piA)^ ^ Vt T / . (5) 

TT £^ y4(d - 1) - '\2y^{d^)J 

The parameters yl-^'' are defined as 

with y}^^ the number of t-pcriodic walks where no back scattering is allowed (nb walks). 
Since (^Y^"^^^ = [d — 1)*, y['^^ is the properly regularized deviation of the number of t- 
periodic nb-walks from their mean. In combinatorial graph theory it is customary to 

define Ct = as the number of t-periodic nb cycles. It is known that for t < \og^_i V 
and asymptotically in V, the Ct are distributed like independent Poisson variables 
[6, 7, 8, 9]. 
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The trace formula for the spectral density of magnetic graphs is similar to (5) 
with the following differences: The entire spectrum of the magnetic adjacency matrix is 
included in the definition of the spectral density. The coefficients y^.^^ are now defined 
as 

(M) _ i- Yt y{M) _ ix„ (r,^ 

where the sum above is over all the nb t-periodic walks, and Xa is the total phase (net 
magnetic flux) accumulated along the t-periodic walk. For finite V , (yi^^^ 7^ since 
Xa = for nb periodic walks where each bond is traversed equal number of times in the 
two directions. However, in the limit of large V, the number of such walks is small and 
therefore (^[^"^^ 0. 

1.2. Spectral fluctuations on graphs and RMT- Numerical evidence 

So far, the only evidence suggesting a connection between RMT and the spectral 
statistics on graphs is the numerical studies of Jacobson et. al. [2]. In a preliminary 
step in the present research, we performed numerical simulations which extended the 
tests of [2] . While describing these studies, we shall introduce a few concepts from RMT 
which will be used in the main body of the paper. 

It is advantageous to map the spectrum from the real line to the unit circle, 

4>j = arccos — < 4>j < ir . (8) 
2\'d — 1 

This change of variables is allowed since in the limit of large graphs, only a fraction 
of order 1/V of the spectrum is outside the support of the Kesten McKay distribution 

[-2y/d^,2Vd^] [10]. 

The mean spectral density on the circle is not uniform, and the Kesten McKay 
density on the circle is 

. U\ 2((i-l) sin^0 

P.M(0) = ^^^-^(py^. (9) 
The mean spectral counting function is defined as 

Nkm{(I>) pKM{(t>)d(t> = - '^-j^ arctan (^^^3^ tan 0^ ^ . (10) 

Following the standard methods of spectral statistics, one introduces a new variable 
9^ which is uniformly distributed on the unit circle. This "unfolding" procedure is 
explicitly given by 

= yNKM{<t>i) (11) 

The nearest spacing distribution defined as 
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is often used to test the agreement with the predictions of RMT (This was also the test 
conducted in ([2])). In this definition of the nearest spacing distribution, 6^0 coincides 
with 6v, since the phases he on the unit circle. In figure (1) we show numerical 
simulations obtained by averaging over 1000 randomly generated 3-regular graphs on 
1000 vertices and their "magnetic" counterparts, together with the predictions of RMT 
for the COE and the CUE ensembles [3], respectively. The agreement is quite impressive. 



in 





Figure 1. Nearest level spacings for: (a.) Graphs possessing time reversal symmetry, 
(b.) Magnetic graphs. 

Both figures are accompanied with the RMT predictions: Solid line - COE, Dashed 
line - CUE. 



Another quantity which is often used for the same purpose is the spectral form- 
factor, 

Kvit) = ^ (^Y.e^'<^' Y ^^^^ 

The form-factor is the Fourier transform of the spectral two point correlation function 
and it plays a very important role in the understanding of the relation between RMT 
and the quantum spectra of classically chaotic systems [3, 13]. 



In RMT the form factor displays scaling: lim Kv{t) = K{t 
limiting expressions for the COE and CUE ensembles are [3]: 

2r - rlog(2r + 1), for r<l 



V 



Kcoe{t) 



). The explicit 



(14) 



rlog 



2t+1 
2t-1' 



for r > 1 



Trace Formulae and Spectral Statistics for Discrete Laplacians on Regular Graphs (11)5 



{r, for r < 1 
(15) 
1, for r > 1 

The numerical data used to compute the nearest neighbor spacing distribution 
P(s), was used to calculate the corresponding form factors for the non- magnetic and the 
magnetic graphs, as shown in figure (2). The agreement between the numerical results 
and the RMT predictions is apparent. This numerical data triggered the research which 
is reported in the present article. 





Figure 2. The form factor K{t) (unfolded spectrum) for (a.) 3-regular graphs 
numerical vs. the COE prediction, (b.) 3-regular magnetic graphs vs. the CUE 
prediction. 



The above comparisons between the predictions of RMT and the spectral statistics 
of the eigenvalues of ci-regular graphs was based on the unfolding of the phases (pj into 
the uniformly distributed phases 6j. As will become clear in the next sections, it is 
more natural to study here the fluctuations in the original spectrum and in particular 
the form factor 

1 ' 



Kv{t) 



V 



The transformation between the two spectra was effected by (11) which is one-to-one 
and its inverse is defined: 



<p = s{e) = N-\Av 



2tx 



(17) 
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This relationship enables us to express Ky{t) in terms of Kv{t). In particular, if Ky{t) 
scales by introducing t = ^ then, 

K{r = ^) = ^l\eK{TS'{e)). (18) 

The derivation of this identity is straightforward, and is described in Appendix A. 

Figure (3) shows K{t = y) = Kyit) computed by assuming that its unfolded 
analogue takes the RMT form (14) or (15), and it is compared with the numerical 
data for graphs with d = 10. It is not a surprise that this way of comparing between 
the predictions of RMT and the data, shows the same agreement as the one observed 
previously. 




T 7 

Figure 3. The form factor K{t) (original spectrum) for 10-regular graphs. The 
numerical results are presented vs. the expression (18) assuming RMT in the dashed 
line, and K{t) (14, 15) in the solid line. 

Remark: The rather unusual definition of the form factor for the original spectrum 
can be illustrated by the following example. Consider the Gaussian ensemble of 
N dimensional symmetric matrices (GOE). Its spectrum (properly normalized) is 
supported on the interval [—1,1] and the mean spectral density is given by Wigner's 
semi-circle law. Mapping the spectrum onto the unit circle results in points which 
are non-uniformly distributed. One can generate the form factors K(t) and K{t) 
from the original and unfolded spectra, and compare the numerical distributions to 
the predictions from GOE. The corresponding K{t) is obtained from (18) in the limit 
d > 1. 

With this summary of definitions and numerical data we prepared the background 
for the main results of the present work, where we use the trace formulae to express 



Trace Formulae and Spectral Statistics for Discrete Laplacians on Regular Graphs (11)7 



the spectral form factor in terms of the variance of the fluctuations in the counting of 
t-periodic nb walks on the two graph ensembles. In Chapter 2 we shall use the known 
properties of t-periodic nb walks to compute the leading term in the Taylor series of 
K{t), near r = 0. In Chapter 3 we shall take the opposite direction, and by assuming 
that the spectral fluctuations for the graphs are given by RMT, we shall derive new 
expressions for the counting statistics of t-periodic nb walks on graphs. This approach 
is similar in spirit to the work of Keating and Snaith [12] who computed the mean 
moments of the Riemann C function on the critical line, assuming that the fluctuations 
of the Riemann zeros follow the predictions of RMT for the CUE ensemble. 



2. From counting statistics of ^-periodic orbits to RMT 

In this section we shall establish a rigorous connection between the spectral properties 
of regular graphs and those predicted by RMT. To achieve this goal, we use the trace 
formula (5), where yt are deflned by (6) or (7) for the two ensembles. 

Defining = arccos ^ 2^^ ) ' using the orthogonality of the cosine, we can extract 
yt, 

yt^2 cos {t(f))p{<j))d(f) (20) 

And so: 

{yl)^A / cos{t(t>)cos{ti;){p{(t>)pmd(t>di; . (21) 
Jo Jo 

RecaUing (A.2) 

PIT PTT 

Kv{s) = 2V / cos{s(t))cos{sij){p{(t))p{ij))d4)di} , (22) 
Jo Jo 

and comparing (21) and (22) we get: 

Kv{t) = \ {yl) . (23) 

So far the treatment of the two ensembles was carried on the same formal footing. We 
shall now address each ensemble separately. 

2.1. The form factor for the Qv4 ensemble 

As was mentioned previously, it is useful to define the number of nb t-cycles on the 



graph as Ct = (in this definition, one does not distinguish between cycles which 
are conjugate to each other by time reversal). From combinatorial graph theory it is 
known that on average (Ct) = [6, 7, 8, 9]. Hence K can also be written as: 

Kif\t) = - • ' . (24) 



2t 

2\ 



t_ {{a - { Ct)Y) 
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The expression of the spectral form factor K^\t) in terms of combinatorial quantities 
is the main result of the present work. In particular, it shows that the form factor is the 
ratio between the variance of C^-the number of nb t-cycles - and its mean. This relation 
is valid for all t in the limit V oo. 

For t satisfying t < loga-iV, it is known that asymptotically, for large V, the Cj's 
are distributed as independent Poisson variables. For a Poisson variable, the variance 
and mean are equal. This implies that for r = ^ — >■ 0, Klf^{t) — r. Notice that the 
relation (18) imphes Kv{t) = 2Kv{t) for r ^ (see also A.9). 

Thus, for T < 1 



This result coincides with the COE prediction (14). It provides the first rigorous 
support of the connection, estabhshed so far numerically, between RMT and the spectral 
statistics of graphs . It is analogous to Berry's "diagonal approximation" [13] in quantum 
chaos. 



In the magnetic ensemble the matrices (2) are complex valued and Hermitian, which is 
tantamount to breaking time reversal symmetry. The relevant RMT ensemble in this 
case is the CUE. 

In the ensuing derivation we shall take advantage of the statistical independence 
assumed for the magnetic phases which arc uniformly distributed on the circle. Ensemble 
averaging will imply averaging over both the magnetic phases and the graphs. 

Recall that yt^^\ in the case of magnetic graphs, was defined, by (7): 



is the sum of interfering phase factors contributed by the individual nb ^-periodic 
walks on the graph. The phase factors of periodic walks which are related by time 
reversal arc complex conjugated. Periodic walks which arc self tracing (meaning that 
every bond on the cycle is traversed the same number of times in both directions), have 
no phase: = 0. Using standard arguments from combinatorial graph theory one can 
show that for t < log^^i V, self tracing nb t-periodic walks arc rare. Moreover, the 
number of t-periodic walks which are repetitions of shorter periodic walks can also be 
neglected. Hence 



primitive cycles. The number of t-cycles on the graph is Ct, hence (26) has approximately 
Ct terms. From (26) and the definition of yl'^''\ it is easily seen that {{y\'^^Y) = 



(25) 



2.2. The form factor for the Qy^^ ensemble 
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^2(J_]^)t 4t^(^^ cos(xq)^ )■ Averaging over the independent magnetic phases we get 
that 



t T 



For T — > 0, and to leading order, K — 2K (A. 9). Hence 
which agrees with the CUE prediction. 



(27) 



(28) 



3. Counting statistics of ^-periodic cycles on d-reguleir graphs from RMT 

In the previous section we made use of the known asymptotic statistics of Yf to show 
that the leading term in the expansion of Kv{t) behaves as gr where = 1, 2 for the 
two graph ensembles. This property is consistent with the predictions of RMT. Had 
we known more about the counting statistics, we could make further predictions and 
compare them to RMT results. However, to the best of our knowledge we have exhausted 
what is known from combinatorial graph theory, and the only way to proceed would be 
to take the reverse approach, and assume that the form factor for graphs is given by 
the predictions of RMT, and see what this implies for the counting statistics. Checking 
these predictions from the combinatorial point of view is beyond our scope. However, 
we shall show that they are accurately supported by the numerical simulations. 

The starting points for the discussion are the relations (23) and (18) which can be 
combined to give 



{{ytf)-lKy{t)^^jyeK(rS'{9)) . 



(29) 



V ''^ ' Vt: 

Our strategy here will be to use for the unfolded form factor the known expressions 
from RMT (14) and (15) and compute ((y*)^)- This will provide an expression for the 
combinatorial quantities defined for each of the graph ensembles, and expanding in r 
we shall compute the leading correction to their known asymptotic values. 

To proceed, we have to analyze the integral (18) and expand it in powers of r near 
r = 0. For this purpose we have to recall the function S{6) = N^\j (^^) i^"^)- The 
inversion of the spectral counting function needs more attention near the end points of 
the support, where 



2V 
<t>-^o 67r 



D-= 

Thus, in the vicinity of ^ = 



S{9) ^ So{9) 



d{d-l) 



30 
AD 



and S' which is singular near is 



S'{9) ^ S',{9) = 



AD 



AD 
'39 



(30) 



(31) 



(32) 
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Both RMT form factors take the value 1 for argument values sufficiently larger than 1. 
Hence the value of 6 where the tS'{9) = 1 plays an important role. We shall denote it 
by Om{T), and for small values of r it takes the value 



9m{T) = 



T 



ADi 



This information suffices for the analytic derivations which will follow. However, 

for the numerical computation of K{t) for the entire range of r, we shall need a better 
approximation for S{6) valid over the entire range of integration. This can be achieved 
by successive Newton-Raphson iterations. To second order, 

3(6) = -^o+^^p^ ^26 - 2dSo + 2{d - 2) arctan . (^^2 ^ _ ^^^2 ^^^^^ ^33) 

The numerical error, induced by this approximation, is less than 5 percents, (it is less 
than 2 percents for > 8). 

We now turn to the small r domain. We shall treat the two ensembles separately 
starting with the magnetic ensemble since it is simpler. 

3.1. Counting statistics for the Qy^^ ensemble 

The CUE form factor (15) is: Kcue{t) = t for r < 1, and Kcue{t) = 1 for 1 < t. 
Therefore, we can divide the integral (18), in the following way: 

^(M)(^) = \ {^fj do + j^^ TS\e)de^ ^1(0^ + "^- TS{em)) =(34) 

= J + - (^m - TS{em)) 
Z TT 

Hence, the first two terms in the expansion of K{t) are 



K^''\r) = ''- + fi{d)Ti + ... (35) 



where 



Thus, the difference (K^^\t) — |)//i((i), should scale for small r as for all values 
of d. This data collapse is shown in figure (4). 

3.2. Counting statistics for the Qv,d ensemble 

Similar results can be obtained for the counting statistics of the Qy^^ ensemble. Here, 
the relevant RMT form factor is the COE expression (14). The integral (18), is divided 
in the following way: 
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Figure 4. (X(*^)(r) - |)//i(d) for various values of d vs. the curve . 



1 

TT 



2tS'{9) - tS'{9) log {2rS'{9) + l)dd 



26^ ^TiT- 2TS{em) - r / 5' [e) log ^ . ci^ 



TT 



2r^'(^) - 1 



7r 



+ - -r / 5>) log (2r5>) + \)de 



(38) 
(39) 
(40) 



We change variables to integrate over the new variable x = S{9). Denoting S' (9) by 
f{x), We get: 



1 



1 / fSiem) 2rf(r) + 1 

20„, + vrr - 2r5(0^) - r / rfx log ^ ^ 1 (41) 



TT 



2r/(x) - 1 



+ - 



TT 



r / cia;log(2r/(x) + 1) 

S(6»„) 



At small values of r (which imply small values of 6m), f{x) = The first integral 

in (41) can be solved explicitly. In the second integral, we can restrict ourselves only to 
the interval in which f{x) oc ^ because the rest of the integral will only yield higher 
order terms in r. Therefore, we can also solve the second integral. Finally we get 

i^(-4)(r) = r + /2(rf)ri + ... (42) 

where 

1 



f2{d) 



'2D \ TT 



■ arccoth(v^) — 1 



(43) 
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Thus, the difference (K'^^^t) — t)/ f2{d), at small r, should scale as independently 
of d. This data collapse is shown in figure (5). 




Figure 5. (X(^)(t) -r)//2(d) for various values of d vs. the curve r 2 . 



Using the above and (24) we can write 

TIT 



dOK, 



COE 



rs\e))-~^i + Ud)V^ + ... 

J T— !>U 



(44) 



If the Ct's were Poissonian random variables, the expansion above would terminate at 
1. Since it does not, we must conclude that the Ct's are not Poissonian. The highest 
order deviation comes from the next order term in the expansion which is proportional 
to r5. The coefficient, /2((i), is explicitly calculated above. 

We can examine the behavior at another domain of r, namely r > 1. 
It can easily be shown that S' {Q) > Consequently, for r > the argument 

of K in (18) is larger than one, and so 

lim ^(^)(r) = 1. (45) 

T— >00 

Combining this result with (24), provides the asymptotic of the variance-to-mean ratio: 
_var{Ct) 



lim r- 



1, for V, t — )• 00; 



± 

V 



T 



(46) 



This is a new interesting combinatorial result, since very little is known about the 
counting statistics of periodic orbits in the regime of r > 1. 



Trace Formulae and Spectral Statistics for Discrete Laplacians on Regular Graphs (11)13 
4. Discussion 

The two main results of the present paper can be summarized as follows. First, we 
have shown that to leading order, the spectral statistics for graphs with time reversal 
symmetry is consistent with the COE, and when time reversal symmetry is broken, the 
CUE statistics come to play. 

Second, by inverting the argument, and assuming RMT for the spectral statistics, we 
derived new results in graph theory, namely the deviation of the number of cycles from 
complete randomness, and the statistics at large r . We do not know at this point how 
to interpret these results from a combinatorial point of view. This remains for now 
an open question. It is important to emphasize that unlike the standard approach in 
RMT, in this paper we have worked with the entire spectrum (bulk and edge states), 
not merely the bulk. As a result, effects of the edge states must be taken into account 
when trying to give a combinatorial answer to the questions posed above. 
So far our rigorous results are rather limited. Yet, this work paves the way to further 
studies where the intricate relationship between combinatorial graph theory and RMT 
will be elucidated. 
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Appendix A. The relation between K{t) and i^(T) 
The form factor (16) can be rewritten in the form 



Jo Jo 

The factor V above is due to the normahzation (1) of the spectral density. The smooth 
part of the spectral density does not encode any information about spectral fluctuations. 




(A.l) 
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so we are only interested in the fluctuating part, p. Thus, 

Kv{s) = 2V / cos{s<f))cos{stlj){p{<j))p{ip))d(i)d^ (A.2) 
Jo Jo 

We emphasize again that the main difference between K and the actual form factor, 
comes from the fact that the 0's are not uniformly distributed. Using the mapping 
(l)^S{e) (17), we get: 

Kvit)=2V / cos{tS{9))cos{tS{9')){p{9)p{e'))d9d9' (A.3) 
Jo Jo 

The two-point correlation function is defined as R2{w) — {2t:)^ (^p{9)p{9 -\- w"^)) . 
In addition, we change variables to 77 = w — 9 — 9' , and we expand the integrand 
up to first order in w, keeping in mind that R2{w) is only of significant magnitude if w 
is small. We are thus left with: 

Kv{t) = ^ j^J dr] dwR2 (^-^^ ■ [cos{2tS{r])) + cos{twS' {r]))j (A.4) 

The first integral is: 

— I drjj ^ dwR2 i^—j cos{2tS{v)) (A.5) 

-1 /■27r POO 

= — / dr) COS {2tS{r))) dsR2is) ^ S{t) (A.6) 
27r Jo J -00 

where we used f^^dsR2{s) = 1 (see for example ([11])). 
The second integral is: 



V 



4^2 dwR2 (^1^) cos{twS'{rj)) = ^ d^K (tS'{v)) (A.7) 

And we conclude that: 



Kv{t) = 5{r) + l-J^\nK (rS'iri)) (A.8) 

Where r is defined as before, and Kv{t) admits the same scaling as in RMT: K{r = 
t/V). Finally, we drop the 5-function term, and we take advantage of the fact that 
5" is symmetric around tt (this is a consequence of the Kesten-McKay measure being 
symmetric around zero). This completes the proof, and we end up with (18): 



Kv{t) = IJ^ d9K (tS'{9)^ 



Using this relation, we can prove that the slope of K is twice that of X at r = 0. 
Denote = g. Then, 

dk{r = 0)_g rs'{9)d9^'- (A.9) 



dr n Jn 2 

which proves the above. 
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